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ABSTRACT 

We use a Bayesian software package to analyze CARMA-8 data towards 19 unconfirmed 
Planck SZ-cluster candidates from Rodrfguez-Gonzalvez et al. (2015) that are associated 
with significant overdensities in WISE. We used two cluster parameterizations, one based 
on a (fixed shape) generalized-NFW pressure profile and another based on a /j gas den¬ 
sity profile (with varying shape parameters) to obtain parameter estimates for the nine 
CARMA-8 SZ-detected clusters. We find our sample is comprised of massive, (Isoo) = 
0.0010 ± 0.0015 arcmin 2 , relatively compact, ( 6500 ) = 3.9 ± 2.0' systems. Results from the /3 
model show that our cluster candidates exhibit a heterogeneous set of brightness-temperature 
profiles. Comparison of Planck and CARMA-8 measurements showed good agreement in 
F 500 and an absence of obvious biases. We estimated the total cluster mass M 500 as a function 
of z for one of the systems; at the preferred photometric redshift of 0.5, the derived mass, 
M 500 ~ 0.8 ± 0.2 x 1O 15 M 0 . Spectroscopic Keck/MOSFIRE data confirmed a galaxy mem¬ 
ber of one of our cluster candidates to be at z = 0.565. Applying a Planck prior in F 500 to 
the CARMA-8 results reduces uncertainties for both parameters by a factor > 4, relative to 
the independent Planck or CARMA-8 measurements. We here demonstrate a powerful tech¬ 
nique to find massive clusters at intermediate (z > 0.5) redshifts using a cross-correlation 
between Planck and WISE data, with high-resolution follow-up with CARMA-8. We also use 
the combined capabilities of Planck and CARMA-8 to obtain a dramatic reduction by a factor 
of several, in parameter uncertainties. 


1 INTRODUCTION 

The Planck satellite (Tauber et al. 2010 and Planck Collaboration 
et al. 2011 I) is a third generation space-based mission to study the 
Cosmic Microwave Background (CMB) and its foregrounds. It has 
mapped the entire sky at nine frequencies from 30 to 857 GHz, with 
an angular resolution of 33 r to 5', respectively. Massive clusters 
have been detected in the Planck data via the Sunyaev-Zel’dovich 
(SZ) effect (Sunyaev & Zel’dovich 1972). Planck has published a 
cluster catalog containing 1227 entries, out of which 861 are con¬ 
firmed associations with clusters. 178 of these were previously un¬ 
known clusters while a further 366 remain unconfirmed (Planck 
Collaboration et al. 2013 XXIX). The number of cluster candidates 
identified in the second data releases, PR2, has now reached 1653. 
Eventually, the cluster counts will be used to measure the cluster 
mass function and constrain cosmological parameters (Planck Col¬ 
laboration et al. 2013 XX). However, using cluster counts to con¬ 


strain cosmology relies, amongst other things, on understanding the 
completeness of the survey and measuring both the cluster masses 
and redshifts accurately (for a comprehensive review see, e.g., Voit 
2005 and Allen, Evrard, & Mantz 2011). To do so, it is crucial to 
identify sources of bias and to minimize uncertainty in the transla¬ 
tion from cluster observable to mass. Regarding cluster mass, since 
it is not a direct observable, the best mass-observable relations need 
to be characterized in order to translate the Planck SZ signal into a 
cluster mass. 

The accuracy of the Planck measurements of the integrated SZ 
effect at intermediate redshifts where, e.g., X-ray data commonly 
reach out to, is limited by its resolution (~ 10' at SZ-relevant fre¬ 
quencies) because the integrated SZ signal exhibits a well-known 
degeneracy with the cluster angular extent (see e.g., Planck Col¬ 
laboration et al. 2013 XXIX). Higher resolution SZ follow-up of 
Planck -detected clusters can help constrain the cluster size by mea¬ 
suring the spatial profile of the temperature decrement and identify 
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sources of bias. Moreover, a recent comparison of the integrated Throughout this work we use J2000 coordinates, as well as a 

SZ signal measured by the Arcminute MicroKelvin Imager (AMI; ACDM cosmology with £2 m = 0.3, £2 a = 0-7, = 0, flj, = 0.041, 

AMI Consortium: Zwart et al. 2008) on arcminute scales and by w 0 = — 1, w a = 0 and tfg = 0.8. Hq is taken as 70kms~ 1 Mpc . 

Planck showed that the Planck measurements were systematically 
higher by « 35% (Planck Collaboration et al. 2013e). This study, 
and its follow-up paper on 99 clusters (AMI and Planck Consor¬ 
tia: Perrott et al. 2014), together with another one by Muchovej et 
al. (2012) comparing CARMA -8 and Planck data towards two sys¬ 
tems, have demonstrated that cluster parameter uncertainties can be 
greatly reduced by combining both datasets. 

In this work, we have used the eight element C ARMA interfer¬ 
ometer, CARMA -8 (see Muchovej et al. 2007 for further details) to 
undertake high spatial resolution follow-up observations at 31 GHz 
towards 19 unconfirmed Planck cluster candidates 1 (Table 1). Our 
primary goal was to attempt to identify massive clusters at high 
redshifts. For this reason, our candidate clusters were those Planck 
SZ-candidates that had significant overdensities of galaxies in the 
WISE early data release (Wright et al. 2010) (>1 galaxy/arcmin 2 ) 
and a red object 2 within 2.5' fainter than 15.8 Vega magnitudes in 
the WISE 3.4-micron band (which corresponds to a 10L* galaxy at 
z ~ l 3 ), to maximize the chances of choosing z > 1 systems. Sim¬ 
ilar work using WISE to find distant clusters has been undertaken 
by the Massive and Distant Clusters of WISE Survey (MaDCoWS), 
which in Gettings et al. (2012) confirmed their first z ~ 1 cluster. 

This work is presented as a series of two articles. The first 
one, Rodriguez-Gonzalvez et al. (2015), henceforth Paper 1, fo¬ 
cused on the sample selection, data reduction, validation using an¬ 
cillary data and photometric-redshift estimation. This second paper 
is organized as follows. In Section 2 we describe the cluster pa- 
rameterizations for the analysis of the CARMA -8 data and present 
cluster parameter constraints for each model. In addition, we in¬ 
clude Bayesian Evidence values between a model with a cluster 
signal and a model without a cluster signal to assess the quality of 
the detection and identify systems likely to be spurious. Planck- 
derived cluster parameters and estimates of the amount of radio¬ 
source contamination to the Planck signal are given in Section 3. 

Improved constraints in the T 500 — 6500 plane from the application 
of a Planck prior on F 500 to the CARMA -8 results are provided in 
Section 3.3. In Section 4 we discuss the properties of the ensem¬ 
ble of cluster candidates, including their location, morphology and 
cluster-mass estimates and present spectroscopic confirmation for 
one of our targets. In this section we also compare the Planck and 
CARMA -8 data and show how our results relate to similar studies. 

We note that, for homogeneity, since not all the cluster candidates 
in this work are included in the PSZ (Union catalog; Planck Col¬ 
laboration 2013 XXIX), we assign a shorthand cluster ID to each 
system (see Table 1). 


1 The Planck SZ catalog used for the initial selection was an intermediate 
Planck data product known internally as DX7. Planck data are collected and 
reduced in blocks of time. The DX7 maps used in this analysis correspond 
to the reduction of Planck data collected from 12 th August 2009 to the 28 th 
of November 2010, which is the equivalent to 3 full all-sky surveys, using 
the v4.1 processing pipeline. The DX7 maps used in this work are part of 
an internal release amongst the Planck Collaboration members and, thus, is 
not a publicly available data product. It should be noted that the public DR1 
PSZ catalog supersedes the preliminary DX7 catalog used for our selection. 

2 We describe as red , objects whose [3.4]-[4.6] WISE colors are > —0.1 
(in AB mags, 0.5 in Vega). This is known as the Mid-InfraRed (MIR) cri¬ 
terion and has been shown by e.g., Papovich (2008) to preferentially select 
Z > 1 objects. 

3 WISE is sensitive to a galaxy mass of 5 x 10 11 M© at z I. 
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Table 1 . Details of the sample of Planck -detected cluster candidates analyzed in this work and their CARMA-8 observations; for additional information, the 
reader is strongly encouraged to look at Paper 1. For simplicity and homogeneity in the cluster naming convention we use a shorthand ID for the targets. The 
PSZ (Union catalog) name (Planck Collaboration 2013 XXIX) is provided, where available*. The Right Ascension (RA) and Declination (Dec) coordinates 
(J2000) correspond to the map centre of our observations while the cluster coordinates are in Table 4. For the short and long baseline data we provided the 
visibility noise. Targets that have been detected in the CARMA-8 data have their ID highlighted. The SZ decrement towards POM had a low signal-to-noise- 
ratio (SNR) (4.2) and the SZ signal in the CARMA-8 imaged data was considered tentative in Paper 1. 


Cluster 

ID 

Union Name 

RA 

hh mm ss 

Dec 

deg min sec 

Short Baseline (0-2kA) 

0 

(mJy/beam) b 

Long Baseline (2-8 kA) 
o 

(mJy/beam) b 

P014 

PSZ1G014.13+38.38 

16 03 21.62 

03 19 12.00 

0.309 

0.324 

P028 

PSZ1G028.66+50.16 

15 40 10.15 

17 54 25.14 

0.433 

0.451 

P031 

- 

15 27 37.83 

20 40 44.28 

0.727 

0.633 

P049 

- 

14 44 21.61 

31 14 59.88 

0.557 

0.572 

P052 

- 

21 19 02.42 

00 33 00.00 

0.368 

0.386 

P057 

PSZ1G057.71+51.56 

15 48 34.13 

36 07 53.86 

0.451 

0.482 

P086 

PSZ1G086.93+53.18 

15 13 53.36 

52 46 41.56 

0.622 

0.599 

P090 

PSZ 1G090.82+44.13 

16 03 43.65 

59 11 59.61 

0.389 

0.427 

P097 

- 

14 55 13.99 

58 51 42.44 

0.653 

0.660 

P109 

PSZ1G109.88+27.94 

18 23 00.19 

78 21 52.19 

0.562 

0.517 

P121 

PSZ1G121.15+49.64 

13 03 26.20 

67 25 46.70 

0.824 

0.681 

PI 34 

PSZ1G134.59+53.41 

11 51 21.62 

62 21 00.18 

0.590 

0.592 

P138 

PSZ1G138.11+42.03 

10 27 59.07 

70 35 19.51 

2.170 

0.982 

P170 

PSZ1G171.01+39.44 

08 51 05.10 

48 30 18.14 

0.422 

0.469 

P187 

PSZ1G187.53+21.92 

07 32 18.01 

31 38 39.03 

0.411 

0.412 

P190 

PSZ1G190.68+66.46 

11 06 04.09 

33 33 45.23 

0.450 

0.356 

P205 

PSZ1G205.85+73.77 

11 38 13.47 

27 55 05.62 

0.385 

0.431 

P264 

- 

10 44 48.19 

-17 31 53.90 

0.476 

0.513 

P351 

- 

15 04 04.90 

-06 07 15.25 

0.355 

0.392 


(t) Since the cluster selection criteria, as well as the data for the cluster extraction, are different to those for the PSZ catalog, not all the clusters in this work have an official Planck ID. 
(a) Achieved rms noise in corresponding maps. 


2 QUANTITATIVE ANALYSIS OF CARMA-8 DATA 
2.1 Parameter Estimation using Interferometric Data 


through the likelihood function (see Feroz et al. 2009 for a detailed 
overview). 

The observed SZ surface brightness towards the cluster elec¬ 
tron reservoir can be expressed as 


In this work we have used McAdam, a Bayesian analysis pack¬ 
age, for the quantitative analysis of the cluster parameters. This 
package has been used extensively to analyze cluster signals in in¬ 
terferometric data from AMI (see e.g., AMI Consortium: Scham- 
mel et al. 2012, AMI Consortium: Rodriguez-Gonzalvez et al. 
2012 & AMI Consortium: Shimwell et al. 2013 for real data and 
AMI Consortium: Olamaie et al. 2012 for simulated data) and 
once before on CARMA-8 data (AMI Consortium: Shimwell et 
al. 2013b). McAdam was originally developed by Marshall, Hob¬ 
son, & Slosar (2003) and later adapted by Feroz et al. (2009) 
to work on interferometric SZ data using an inference engine, 
MultiNest (Feroz & Hobson 2008 & Feroz, Hobson & Bridges 
2008), that has been optimized to sample efficiently from com¬ 
plex, degenerate, multi-peaked posterior distributions. McAdam 
allows for the cluster and radio source/s (where present) param¬ 
eters to be be fitted simultaneously directly to the short baseline 
(SB; ~ 0.4 —2 kA) uv data in the presence of receiver noise and 
primary CMB anisotropies. The high resolution, long baseline (LB; 
~ 2 — 10 kA) data are used to constrain the flux and position of de¬ 
tected radio sources; these source-parameter estimates are then set 
as priors in the analysis of the SB data (see Section 2.2.3). Our 
short integration times required all of the LB data to be used for the 
determination of radio-source priors and none of the LB data were 
included in the McAdam analysis of the SB data. Undertaking the 
analysis in the Fourier plane avoids the complications associated 
with going from the sampled visibility plane to the image plane. In 
McAdam, predicted visibilities Vy (u,) at frequency v and base¬ 
line vector u ( , are generated and compared to the observed data 


A/cmb = A7cmb 


dB(v.T) 

dT 


Tcmb 


0 ) 


where —yrp |r CM B derivative of the black body function at 

Tcmb - the temperature of the CMB radiation (Fixsen et al. 1996). 
The CMB brightness temperature from the SZ effect is given by 


A7cmb = /(v)yTcMB- (2) 

Here, /(v) is the frequency (v)-dependent term of the SZ effect, 

m={f^\-^(i + Ssz(x,T e )), (3) 

where the §sz term accounts for relativistic corrections (see Itoh 
et al. 1998), T e is the electron temperature, x = hv /kg7cMB, h is 
Planck’s constant and kg is the Boltzmann constant. To calculate 
the contribution of the cluster SZ signal to the (predicted) visibil¬ 
ity data, the Comptonization parameter, y, across the sky must be 
computed: 

yW = — «e(r)k B T e (r)dloc f P e (r)dl. (4) 

Here, Or is the Thomson scattering cross-section, m e is the electron 
mass, n e (r), T e (r) and P e (r) are the electron density, temperature 
and pressure at radius r respectively, c is the speed of light and dl is 
the line element along the line of sight. The projected distance from 
the cluster center to the sky is denoted by s, such that i~ = s 2 + l 2 . 
The integral of y over the solid angle dfl subtended by the cluster 
is proportional to the volume-integrated gas pressure, meaning this 
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Table 2. Summary of the cluster priors used in our analysis for the observa¬ 
tional gNFW cluster parameterization (model I, Section 2.2.1). Ax c and Ay c 
are the displacement from the map center to the centroid of the SZ decre¬ 
ment in RA and Dec, respectively, rj is the ellipticity parameter and the 
position angle. 6 S — r s /D\, where r s is the scale radius and D\ is the angu¬ 
lar size distance to the cluster. Yi is the SZ surface brightness integrated over 
the cluster’s extent on the sky. The 6 S and Yj priors have been previously 
used in Planck Collaboration et al. (2011 VII) and Planck Collaboration et 
al. (2013e). 


Parameter 

Prior 

Ax c 

Gaussian centered at pointing centre, <7 = 60 // 

Ay c 

Gaussian centered at pointing centre, <7 = 60" 

n 

Uniform from 0.5 to 1.0 

<f> 

Uniform from 0 to 180° 

e s 

^_g-A0s f or 

1.3' < 0j < 45' & 0 outside this range 

Tt 

Yj 01 for 0.0005 to 0.2arcmin 2 
& 0 outside this range with a — 1.6 


quantity correlates well with the mass of the cluster. For a spherical 
geometry this is given by 

^sphM = —^ [ P e (r')4n/ 2 dr' . (5) 

m e c- Jo 

When r —> °°, Equation 5 can be solved analytically, as shown in 
AMI and Planck Consortia: Perrott et al. (2014), yielding the total 
integrated Compton-v parameter, Ir.phys* which is related to the SZ 
surface brightness integrated over the cluster’s extent on the sky 
through the angular diameter distance to the cluster (D A ) as Yj = 
^T,phys/^A' 


2.2 Models and Parameter Estimates 

Analyses of X-ray or SZ data of the intra-cluster medium (ICM) 
that aim to estimate cluster parameters are usually based on a pa¬ 
rameterized cluster model. Cluster models necessarily assume a ge¬ 
ometry for the SZ signal, typically spherical, and functional forms 
of two linearly-independent thermodynamic cluster quantities such 
as electron temperature and density. These models commonly make 
assumptions such as, the cluster gas is in hydrostatic equilibrium or 
that the temperature or gas fraction throughout the cluster is con¬ 
stant. Consequently, the accuracy and validity of the results will de¬ 
pend on how well the chosen parameterization fits the data and on 
the effects of the model assumptions (see e.g., Plagge et al. 2010, 
Mroczkowski 2011 & AMI Consortium: Rodriguez-Gonzalvez et 
al 2011 for studies exploring model effects in analyses of real data 
and AMI Consortium: Olamaie et al. 2012 and Olamaie, Hob¬ 
son, & Grainge 2013 for similar work on simulated data). In this 
work we present cluster parameters calculated from two different 
models; one is based on a fixed-profile-shape gNFW parameteri¬ 
zation, for which typical marginalised parameter distributions for 
similar interferometric data from AMI have been shown in e.g., 
AMI and Planck Consortia: Perrott et al. (2014), and a second is 
based on the jS profile with variable shape parameters, where typi¬ 
cal marginalised parameter distributions for comparable AMI data 
have been presented in AMI Consortium: Rodriguez-Gonzalvez et 
al. (2012). Comparison of marginalised posteriors for CARMA and 
AMI data in AMI Consortium: Shimwell et al. (2013b) for the /3 
model showed the distributions to be very similar. The clusters pre¬ 
sented here are at modest redshifts and are unlikely to be in hydro¬ 
static equilibrium - adopting two models at least allows the depen¬ 


dency of the cluster parameters on the adopted model to be illus¬ 
trated and a comparison with previous work to be undertaken. 


2.2.1 Cluster model I: observational gNFW parameterization 


For cluster model I, we have used a generalized-NFW (gNFW; 
Navarro, Frenk, & White 1996) pressure profile in the same fash¬ 
ion as in the analysis of Planck data (Planck Collaboration et al. 
2011 VII) to facilitate comparison of cluster parameters. A gNFW 
pressure profile with a fixed set of parameters is believed to be a 
reasonable choice since ( 1 ) numerical simulations show low scat¬ 
ter amongst cluster pressure profiles, with the pressure being one 
of the cluster parameters that suffers least from the effects of non- 
gravitational processes in the ICM out to the cluster outskirts and 
( 2 ) the dark matter potential plays the dominant role in defining the 
distribution of the gas pressure, yielding a (pure) NFW form to the 
profile, which can be modified into a gNFW form to account for the 
effects of ICM processes (see e.g., Vikhlinin et al. 2005 and Nagai, 
Kravtsov, & Vikhlinin 2007). Using a fixed gNFW profile for clus¬ 
ter models has become regular practice (e.g., Atrio-Barandela et al. 
(2008) for WMAP. Mroczkowski et al. (2009) for SZA, Czakon et 
al. (2014) for BOLOCAM and Plagge et al. (2010) for SPT data). 

Assuming a spherical cluster geometry, the form of the gNFW 
pressure profile is the following: 


Pe(r)=Po 



1+1 - 

r s 


( c-b)/a 


( 6 ) 


where Pq is the normalization coefficient of the pressure profile 
and r s is the scale radius, typically expressed in terms of the con¬ 
centration parameter C 50 O = r^oo/ r s- Parameters with a numerical 
subscript 500, like C 50 O, refer to the value of that variable within 
r 500—the radius at which the mean density is 500 times the crit¬ 
ical density at the cluster redshift. The shape of the profile at in¬ 
termediate regions (r ~ r s ), around the cluster outskirts (r » r s ) 
and in the core regions (r « r s ) is governed by three parame¬ 
ters a, b, c, respectively. Together with C 500 , they constitute the 
set of gNFW parameters. Two main sets of gNFW parameters 
have been derived from studies of X-ray observations (inner clus¬ 
ter regions) and simulations (cluster outskirts) (Nagai, Kravtsov, 
& Vikhlinin 2007 and Arnaud et al. 2010). For ease of compar¬ 
ison with the Planck results, as well as with SZ-interferometer 
data e.g., from AMI in Planck Collaboration et al. (2013e), we 
have chosen to use the gNFW parameters derived by Arnaud et 
al.: (c 50 o,a,M = (1.156,1.0620,5.4807,0.3292). 

In our gNFW analysis, we characterize the cluster by the fol¬ 
lowing set of sampling parameters (Table 2): 


= (Ax c ,Ay c , t], c E > , 9 S = r s /D A ,Y T ). 


Here, Ax c ,A y c are the displacement of the cluster decrement from 
the pointing centre, where the cluster right ascension is equal to 
the map center (provided in Table 1), rj is the ellipticity parameter, 
that is, the ratio of the semi-minor and semi-major axes and <3> is 
the position angle of the semi-major axis, measured N through E 
i.e. anti-clockwise. We note that the projected cluster decrement is 
modeled as an ellipse and hence our model is not properly triaxial. 

The priors used in this analysis are given in Table 2; they 
have been used previously for the blind detection of clusters in 
Planck data (Planck Collaboration et al. 2011 VII) and to charac¬ 
terize confirmed and candidate clusters in Planck Collaboration et 
al. (2013e). Cluster parameter estimates and the CARMA best-fit 
positions derived from model I are provided in Tables 3 and 4, re¬ 
spectively. 
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Table 3. Mean and 68 %-confidence uncertainties for McADAM-derived cluster parameters when fitting for an observational gNFW cluster parameterization 
(Model I; Section 2.2.1) for clusters with a significant SZ detection in the CARMA -8 data (Table 12). The cluster ID is a shorthand naming convention adopted 
here and in Paper 1, since not all our targets have an identifier in the PLANCK Union catalog (Planck Collaboration et al. 2013 XXIX). Where available, the 
Union catalog names are given in Table 1. The derived sampling parameters for the gNFW parameterization are presented in columns 2 to 7 and their priors are 
listed in Table 2. F 500 is the integrated SZ surface brightness within $ 500 , where $500 = ^ 500 / D& and v(0) is the central Comptonization parameter, y. Equation 
4. 


Cluster ID 

Av c 

Ay c 


rt 

e s 

Ft 

0500 

I 500 

,y(0) 


n 

n 

deg 


/ 

x 10~ 4 arcmin 2 

arcmin 

x 10~ 4 ai-cmin 2 

xl0~ 4 

P014 

25if< 

-148^ 4 

149+g 

°- 6 ±ai 

4^1 

20+}® 

4 3+1- 2 

Hi 

1 6 +0 ' 4 

1 ■ — 0.4 

P086 

68±|j 

91+ 17 

V1 -17 

8 5tg 

0.811 

3-2 

201 

9 7+1.6 

Dm ' -2.2 

Hi 

2 I +0 - 9 
z - A -o.8 

P097 

70+23 

7 —23 

38+|6 

76^52 

0.811 

3+ 1 

Hi 

3 2 +1,2 
J - z -l .3 

61 

1 9 + 0.9 

1 ■ — 0.8 

P109 

10 -Ji 

75+ 14 

7 —13 

89±I4 

0 . 8 l| 

3-1 

io+l 

3 1 +1 

6 1 

1 9 + 0.8 
iy - 0.8 

P170 

-59±Ii 

10+ 14 


0-211 

2 ±! 

12+1 

7 7+0.9 

Z./_i 2 

61 

7 tj+ 0.8 
z ‘ —0.7 

PI 87 

64+ 14 

0+—14 

—67 +15 

0/ —15 

101 +“ 

0.811 

4^1 

79 + 10 
z:) -18 

4 1 + 1-8 

13 +5 

1 9 +O .6 

1 ' — 0.6 

P190 

59LJ8 

nil! 

95+3® 

0-811 

3-1 

l«±ii 

3 3 + 1 - 4 

q+4 

y -6 

2 0+ 06 
z ’ — 0.6 

P205 

-83'!! 

26+Jg 

79±| 

0 7 +0 - 2 

u - 7 -o.i 

4 1 

31+ 13 

Dl -26 

4 9+2-1 
— 2.2 

HI 5 

1 6 +0 - 4 
iO -0.4 

P351 

-42 ±* 

63tM 

7 i +60 

71 -45 

0 7 +0-2 
u - 7 - 0.2 

5±1 

19+ 7 

1 14 

r a+2.3 
^ —2.3 

10 I 

O 

SO 

1 + 

p p 


Table 4. Cluster J2000 coordinates derived using the gNFW (model I) fits 
to the CARMA -8 data. 


Cluster ID 

RA 

hh:mm:ss 

Dec 

dd:mm:ss 

POM 

16:03:23.29 

03:16:44.00 

P086 

15:14:00.85 

52:48:12.56 

P097 

14:55:24.17 

58:52:20.44 

PI 09 

18:23:03.50 

78:23:07.19 

P170 

08:50:59.16 

48:30:28.14 

P187 

07:32:23.03 

31:37:32.03 

PI 90 

11:06:08.81 

33:33:56.23 

P205 

11:38:07.21 

27:54:39.62 

P351 

15:04:02.09 

-06:06:12.24 


Table 5. Summary of the cluster priors used in our analysis for the obser¬ 
vational p cluster parameterization (model II, Section 2.2.2). Ax c and Ay c 
are the displacement from the map center to the centroid of the SZ decre¬ 
ment in RA and Dec, respectively. 77 is the ellipticity parameter and the 
position angle. The power law index p and the core radius r c are the shape 
parameters of the density profile and A7o is the temperature decrement at 
zero projected radius. 


Parameter 

Prior 

Ay c 

Gaussian centered at pointing centre, <7 = 60 7/ 

Ay c 

Gaussian centered at pointing centre, <7 = 60 7/ 

n 

Uniform from 0.5 to 1.0 

<s> 

Uniform from 0 to 180° 

P 

Uniform from 0.4 to 2.5 

0C 

Uniform from 20 to 500" 

ATo 

Uniform from —3 to —0.01 mK 


2.2.2 Cluster model II: observational P parameterization 

For this cluster parameterization we fit for an elliptical cluster ge¬ 
ometry, as we did for model I, and model the shape of the SZ 
temperature decrement with a /3-like profile (Cavaliere & Fusco- 
Femiano 1978): 

1 — 30 

Ar C MB(0)=A7b^l+^j j , (7) 

where ATo is the brightness temperature decrement at zero pro¬ 
jected radius, while /3 and r c = 6 C x D\ — the power law index and 
the core radius—are the shape parameters that give the density pro¬ 
file a flat top at small and a logarithmic slope of 3/5 at large 
The sampling parameters for the cluster signal are: 

= (Ax c ,Ay c ,T7,0,A7o,/3,0 c ), 

with priors given in Table 5, which allow for the y signal to be 
computed. Cluster parameter estimates derived from model II are 
provided in Table 6. 

It is important to note that, historically, in many SZ analyses 
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Table 6. Mean and 68 %-confidence uncertainties for Me ADAM-derived cluster parameters when fitting for an observational /3 cluster parameterization (Model 
II; see Section 2.2.2) for clusters with an SZ detection in the CARMA -8 data (see Table 12). The priors for these sampling parameters are given in Table 5. 
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Figure 1. Left: radial brightness temperature profiles derived from Equation 7 using the cluster parameter values (A7o, /3 and 0 C ) from fits of model II to the 
CARMA-8 data (Table 6). When the profile has dropped to 3 times the noise in the SB data, the lines turn from solid to dashed. In the legend, the cluster 
candidates have been ordered by decreasing F 500 (Table 3), although it should be noted that for some clusters the difference in >500 is small. Shallower profiles 
with the most negative A7o values and largest 6)500 values should correspond to the largest I 500 from model I. Comparing the volume-integrated brightness 
temperature profile, f T(6)47t6 2 d0 for each cluster from 0 to its #500 (from model I; Table 3) with F 500 (model I), shows reasonable correspondence for most 
systems, with the exception of P351 and P187. Right: best-fit radial brightness temperature profiles for three clusters, P205, P014 and P109, which span a 
range of profile shapes (in thick solid lines, same as the corresponding lines in the left panel). However, here we include the upper and lower 68 %-confidence 
limits of each brightness temperature profile in thin solid lines and highlight the region between these limits with diagonal lines. This shows that, according to 
the parameter fits to CARMA -8 data from model II, the clusters in our sample exhibit a heterogenous set of profiles, distinguishable despite the uncertainty. 
Furthermore, high resolution SZ data or X-ray data may be insensitive to the large signal from the outer parts of the cluster and introduce an additional 
uncertainty in the derivation of Yj, an issue which is addressed through a joint analysis with the Planck data. 


the shape of the /3 profile has been fixed to values obtained from fits 
to higher resolution X-ray data. However, these X-ray results pri¬ 
marily probe the inner regions of the cluster and, thus, can provide 
inadequate best-fit profile shape parameters for SZ data extending 
out to larger-r. A comparative analysis in Czakon et al. (2014) re¬ 


veals systematic differences in cluster parameters derived from SZ 
data using a model-independent method versus X-ray-determined 
cluster profiles. Several studies have now shown that fits to SZ data 
reaching r 5 oo and beyond preferentially yield larger /3 values than 
X-ray data, which tend to yield /3 = 2/3 (see e.g., AMI Consor- 
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tium: Hurley-Walker et al. 2012 & Plagge et al. 2010). Using a 
suitable |3 (and 0 C ) value for the aforementioned typical SZ data 
can yield results comparable to those of a gNFW profile. In our ;6 
parameterization we allow the shape parameters, /3 and 0 C , to be fit 
in Me Adam, since they jointly govern the profile shape. While our 
data cannot constrain either of these variables independently, they 
can constrain their degeneracy. 

Although a large fraction of clusters are well-described by 
the best-fit gNFW parameterizations, some are not, as can be seen 
from e.g., the spread in the gNFW parameter sets from fits to in¬ 
dividual clusters in the REXCESS sample (Arnaud et al. 2010). In 
these cases, modelling the cluster using a fixed (inadequate) set of 
gNFW parameters will return biased, incorrect results, whereas us¬ 
ing a |3 model with varying j 8 and 6 C should provide more reliable 
results. This is shown in Figure 1 of AMI Consortium: Rodrfguez- 
Gonzalvez et al. (2012) where data from AMI for a relaxed and 
a disturbed cluster are analyzed with a j 8 parameterization and five 
gNFW parameterizations, four of which have gNFW sets of param¬ 
eters drawn from the Arnaud et al. REXCESS sample, three from 
individual systems and one from the averaged (Universal) profile, 
and, lastly, one with the average-profile values from an independent 
study by Nagai, Kravtsov, & Vikhlinin (2007). For both clusters, the 
Nagai parameterization lead to a larger Igoo — 6500 degeneracy and 
larger parameter uncertainties than the Arnaud Universal parame¬ 
terization. The mean T 500 and d^oo values obtained from using the 
j 8 . Universal (Arnaud) and Nagai gNFW profiles were consistent 
to within the 95% probability contours, but this was not the case 
for fits using the sets of gNFW parameters obtained from individ¬ 
ual fits to REXCESS clusters, indicating that some clusters do not 
follow a single, averaged profile. Here, comparison of the Bayesian 
evidence values for beta and gNFW-based analyses showed that the 
data could not distinguish between them. Our CARMA data for this 
paper have a similar resolution to the AMI data but, typically, they 
have much poorer SNRs and similarly cannot determine which of 
the two profiles provides a better fit to the data. More recently, Say¬ 
ers et al. (2013) have derived a new set of gNFW parameters from 
45 massive galaxy clusters using Bolocam and Mantz et al. (2014) 
have further shown how the choice of model parameters can have a 
measurable effect on the estimated Y-parameter. 

All sets of gNFW parameters can lead to biases when applied 
to different sets of data. Given that there is no optimally-selected set 
of gNFW parameters to represent CARMA 31-GHz data towards 
massive, medium-to-high redshift clusters (z >0.5), we choose to 
base most of our analysis on the gNFW parameter set from the 
’Universal" profile derived by (Arnaud et al. 2010) as this facili¬ 
tates comparison with the Planck analysis and parallel studies be¬ 
tween Planck and AMI, an interferometer operating at 16 GHz with 
areminute resolution. 

2.2.2.1 Cluster Profiles Using Equation 7 for A7cmb($) (the 
SZ temperature decrement), and the mean values for A7o, [5 and 
0 C derived from model II fits to the CARMA -8 data (Table 6 ), in 
Figure 1, left, we plot the radial brightness temperature profiles 
for our sample of CARMA- 8 -detected candidate clusters. We or¬ 
der them in the legend by decreasing CARMA -8 Tsocn from Ta¬ 
ble 3, although in some cases the differences are small. We would 
expect clusters with the most negative A Tq values, the shallowest 
profiles and the largest 6500 to yield the largest >500 values. While 
there is reasonable correspondence throughout our cluster sample, 
two clusters P351 and PI87 are outliers in this relation. Computing 
f T(Q)4jt6 2 dd for each cluster from 0 to its 6500 - determined from 
model I (Table 3) shows that P351 (P187) has the highest (fifth 


highest) volume-integrated brightness temperature profile but only 
the fifth highest (second highest) Y 500 ■ 

In Figure 1, right, we plot the upper and lower limits of the 
brightness temperature profiles allowed by the profile uncertainties 
for three clusters, P014, P109 and P205, chosen to span a wide 
range of profile shapes. It can be seen that the cluster candidates 
display a range of brightness temperature profiles that can be dif¬ 
ferentiated despite the uncertainties. In Table 7, by computing the 
ratio of the integral of the brightness temperature profile within (a) 
the 100-GHz Planck beam and (b) 0 5O o from Table 3, we quantify 
how concentrated each brightness temperature profile is. The pro¬ 
file concentration factors have a spread of a factor of ~ 2.5 but for 
six clusters they agree within a factor of ~ 1.2. Furthermore, the de¬ 
rived ellipticities shown in Table 6 and 3, which can be constrained 
by the angular resolution of the CARMA-8 data, show significant 
evidence of morphological irregularity suggesting that these clus¬ 
ters may be disturbed and heterogeneous systems. 


2.2.3 Radio-source Model and Parameter Estimates 


Radio sources are often strong contaminants of the SZ decrement 
and their contributions must be included in our cluster analysis. In 
this work, we jointly fit for the cluster, radio source and primary 
CMB signals in the SB data. The treatment of radio sources is the 
same for all cluster models. These sources are parameterized by 
four parameters, 

© s = (RA S , Dec s , (Xsi 531 ), 


where RA S and Dec s are RA and Dec of the radio source, a s is the 
spectral index, derived from the low fractional CARMA-8 band¬ 
width and 531 is the 31-GHz integrated source flux. We adopt the 
S v oc v~“ s convention, where S is flux and v frequency. 

The high resolution LB data were mapped in DlFMAP (Shep¬ 
herd 1997) to check for the presence of radio sources. Radio-point 
sources detected in the LB maps were modeled using the DlFMAP 
task MODELFIT. The results from MODELFIT were primary-beam 
corrected using a FWHM of 660 arcseconds by dividing them by 
the following factor: 


exp 



(8) 


where r is the distance of the source to the pointing centre and 

660 

° 2 x (2 x ln(2)) 0 - 5 ’ 

The primary-beam corrected values were used as priors in the anal¬ 
ysis of the SB data (see Table 8). MODELFIT values are given in 
Paper 1 and McADam-derived values are provided in Table 11. 


2.3 Quantifying the Significance of the CARMA-8 SZ 
Detection or Lack thereof 

Bayesian inference provides a quantitative way of ranking model 
fits to a dataset. Although the term model technically refers to a 
position in parameter space ©, here we refer to two model classes'. 
a model class that allows for a cluster signal to be fit to the data, M \, 
and another, Mq, that does not. The parameterization we have used 
for this analysis has been the gNFW-based model, Model I; for the 
M\ case Model I was run as described in Section 2.2.1 and for the 
Mq case it was run in the same fashion except for the prior on Y 500 , 
which was set to 0, such that no SZ (cluster) signal is included in 
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Table 7. Ratio of the integral of J T (0) cl 0 integrated between 0 and 600" 
(f^FWHM of the 100-GHz Planck beam) and integrated between 0 and #500 
(from Table 3), where the expression for T(Q) is given in Equation 7. This 
ratio is a measure of concentration of the brightness temperature profile. 
The most concentrated profile (for P014) is ~ 2.5 times more concentrated 
than the shallowest profile (for P351). Six of the nine clusters are equally 
concentrated to within a factor of « 1.2. 


Cluster ID 

j60°" r ( 0 ) rf0 

j" 500 T(B)dB 

P014 

17.3 

P086 

17.1 

P097 

17.8 

PI 09 

13.5 

PI 70 

36.0 

PI 87 

20.8 

PI 90 

19.4 

P205 

40.5 

P351 

45.07 


Table 8. Summary of the source priors used in our analysis. Values for 
the position (x s , y s ) and 31-GHz flux (S 31 ) priors were obtained from the 
long baseline CARMA -8 data (see Paper 1 and Section 2.2.3 for further 
details). The error on the source location is the typical error of the LB data. 
To account for the fact that in the combined cluster and source analysis 
uses only the lower-resolution data, we model the integrated 31-GHz source 
flux with a Gaussian. The width of this Gaussian is set to <7 = 20% of the 
source flux. For the spectral index a s , we used a wide prior, encompassing 
reasonable values (see Section 3.2). 


Parameter 

Prior 

RA S , Dec s 

Uniform between ±10" 


from the LB-determined position 

S31 

Gaussian centered at best-fit MODELFIT value 


with a a of 20% 

a s 

Gaussian centered at 0.6 with o — 0.5 


the model. Given the data D, deciding whether M\ or Mq fit the data 
best can be done by computing the ratio: 

Pr(M] |D) _ Pr(D|Mi) Pr(MQ _ Z x Pr(MQ 

Pr(M,)|D) Pr(D|M 0 ) Pr(M 0 ) Z 0 Pr(M 0 )' 

Here, K = is known as the Bayes Factor and is the 

prior ratio, that is, the probability ratio of the two model classes, 
which must be set before any information has been drawn from the 
data being analyzed. Here, we set the prior ratio to unity 4 i.e. we 
assume no a priori knowledge regarding which model class is most 
favorable. The Bayesian Evidence Z is calculated as the integral of 
the likelihood function, J?f = Pr(D | ©, M) , times the prior probabil¬ 
ity distribution n(0) = Pr(0|M), 

©,M)Pr(0|M)rf®0 

, (ID 

= / ^f(©)n(©)d^0 

where @ is the dimensionality of the parameter space. Z represents 
an average of the likelihood over the prior and will therefore favour 
models with high likelihood values throughout the entirety of pa¬ 
rameter space. This satisfies Occams razor, which states that the 
models with compact parameter spaces will have larger evidence 
values than more complex models, unless the latter fit the data sig¬ 
nificantly better i.e., unnecessary complexity in a model will be 
penalized with a lower evidence value. 

The derived Bayes factor is listed in Table 12 along with the 
corresponding classification of whether or not the cluster was con¬ 
sidered to be detected. We find that all the SZ decrements con¬ 
sidered to have high signal-to-noise ratios in Paper I have Bayes 
Factors that indicate the presence of a cluster signature is strongly 
favoured. However, we do find some tension between the Paper I 
MODELFIT and McAdam results for one of the candidate clusters, 
P014. In paper I this candidate cluster was catalogued as tentative 
(see Appendix B of Paper I for more details). The low SNR 5 of 4.2 
for the decrement together with the unusually large displacement 
from the Planck position (~ 159”) suggest this detection is spuri¬ 
ous. The lack of an X-ray signature would support this, unless it 
was a high-redshift cluster or one without a concentrated profile. 
With regards to the source environment, two sources were detected 
in the LB data with a peak 31-GHz flux density of 6.3 and 9.4 mJy, 
a distance of ~ 100" and ~ 600" from the SZ decrement, respec¬ 
tively. The LB data, after subtraction of these radio sources using 
the MODELFIT values, were consistent with noise-like fluctuations, 
indicating the removal of the radio-source flux worked well. The 
NVSS results revealed four other radio sources which, due to their 
location and measured fluxes at 1.4 GHz (as well as their lack of 
detection in the LB data), are unlikely to contaminate the candi¬ 
date cluster. The NVSS results also indicate the radio sources are 
not extended. The strongest support for the presence of an SZ sig¬ 
nature comes from the relatively-high Planck SNR of 4.5 but this 
measurement could suffer from the high contamination from inter¬ 
stellar medium emission which mimics itself as the SZ increment 
at high frequencies and could also result in a large error on its de¬ 
rived position - suggestion of this arises from the strength of the 
100/lm emission which is the highest for our sample. Yet, despite 


Z = Pr(D|M) = J Pr(D| 


4 Prior ratios need not be set to unity see e.g, Jenkins & Peacock (2011). 

5 The SNR for the CARMA SZ detections was calculated in Paper I as the 
ratio of the peak decrement, after correcting for beam attenuation, and the 
RMS of the SB data. 
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these results, the Bayes factor from Table 12 shows that a model 
with a cluster signature is preferred over one without. There are 
potentially quite important differences between the Paper I results 
and the McAdam evidences e.g., the MODELFIT results are based 
on an image and single-value fits without simultaneous fits to other 
parameters, while the McAdam results are derived from fits to the 
wv-plane, taking into account all model parameters, some of which 
are not strongly constrained by the data. Indeed, we find high scat¬ 
ter in the relation between evidence values and MODELFIT-based 
SNR values but they are positively correlated. Validation of this 
candidate cluster will require further data. Given the modest sig¬ 
nificance of the detection of P014 by two different techniques, we 
decided to include P014 in our McAdam analyses. 

3 CONSTRAINTS FROM PLANCK 
3.1 Cluster Parameters 

We used the public Planck PR 1 all-sky maps to derive I 500 and 6 5 qq 
values for our cluster candidates (Table 9). The values were derived 
using a multifrequency matched filter (Melin et al. 2006, 2012). 
The y profile from Equation 4 is integrated over the cluster profile 
and then convolved with the Planck beam at the corresponding fre¬ 
quency; the matched filter leverages only the Planck high frequency 
instrument (HFI) data between 100-857 GHz because it has been 
seen that the large beams at lower frequencies result in dilution of 
the temperature decrement due to the cluster. The beam-integrated, 
frequency-dependent SZ signal is then fit with the scaled matched 
filter profile from Equation 2 to derive Tsoo- The uncertainty in the 
derived T 500 is due to both the uncertainty in the cluster size (Planck 
Collaboration et al. 2013 XXIX), as well as the signal to noise of 
the temperature decrement in the Planck data. The large beam of 
Planck , FWHM~ 10' at 100 GHz, makes it challenging to constrain 
the cluster size unless the clusters are at low redshift and thereby 
significantly extended. For this reason, Planck Collaboration et al. 
(2013 XXIX) provided the full range of T 500 — P 500 contours which 
are consistent with the Planck data. 

For the comparison here, there are two Planck- derived T 500 
estimates; T 500 was calculated using the cluster position and size 
(6500) obtained from the higher resolution CARMA -8 data, while 
T500,blind was computed using the Planck data alone without using 
the CARMA -8 size constraints. Similarly, 6500 ,blind is a measure of 
the angular size of the cluster using exclusively the Planck data; this 
value is weakly constrained and, thus, no cluster-specific errors are 
given for this parameter in Table 9. At this point, it is important to 
note that the quoted uncertainty for Tsoo.blind is an underestimate; 
the quoted error for this parameter is based on the spread in T500 
at the best-fit ©soo.blind and is proportional to the signal to noise 
of the cluster in the Planck data i.e, without considering the error 
on 05 oo blind which is very large. However, the uncertainty in *500 
is accurate since it propagates the true uncertainty in 6590 from the 
CARMA -8 data into the estimation of this quantity from the Planck 
maps. 

The uncertainty in Tsoo.blind front using the CARMA-8 size 
measurement has gone down by ~ 60% on average, despite the fact 
that the Tsoo.blind does not include the uncertainty resulting from 
the unknown cluster size. If the true uncertainty in Tsoo.blind had 
been taken into account, the uncertainty would have gone down by 
more than an order of magnitude after application of the CARMA-8 
derived cluster size constraints. The mean ratio of Tsoo.blind to T 500 
is 1.3 and, in fact, T 500 is only larger than its blind counterpart for 


two systems. Differences in the profile shapes account for Tsoo.blind 
being larger than T 500 for three systems, P014, P097 and PI87, for 
which Psoo.blind is smaller than 6500 measured by CARMA- 8 . 

3.2 Estimation of Radio-source Contamination in the Planck 
143-GHz Data 

In order to assess if there are any cluster-specific offsets in the 
Planck T 500 values, we estimate the percentage of radio-source 
contamination to the Planck SZ decrement at 143 GHz—an impor¬ 
tant Planck frequency band for cluster identification—from the 1.4- 
GHz NVSS catalog of radio sources. Spectral indices between 1.4 
and 31 GHz were calculated in Table 3 in Paper 1 for sources de¬ 
tected in both our CARMA -8 LB data and in NVSS, giving a mean 
value of a s of 0.72. We use this value for a s to predict the source- 
flux densities at 100 and 143 GHz of all NVSS sources within 5'of 
the CARMA -8 pointing center, following the same relation as we 
did earlier, S v °c v~“ s . 

The accuracy of the derived 100 and 143-GHz source fluxes 
is uncertain. Firstly, there is source variability due to the fact the 
NVSS and CARMA -8 data were not taken simultaneously, which 
could affect the 1.4-31 GHz spectral index. Secondly, we assume 
the spectral index between 1.4 and 31 GHz is the same as for 1.4 
to 143 GHz, which need not be true. Thirdly, we deduce a s from 
a small number of sources, all of which must be bright in the LB 
data and apply this a s to lower-flux sources found in the deeper 
NVSS data, for which a s might be different. However, previous 
work shows that this value for a s is not unreasonable. Comparison 
of 31-GHz data with 1.4-GHz data on field sources has been previ¬ 
ously done by Muchovej et al. (2010) and Mason et al. (2009). For 
the former, the 1.4-to-31 GHz spectral-index distribution peaked at 
0.7 while, for the latter, it had a mean value of 0.7. The Muchovej et 
al. study also investigated the spectral index distribution between 5 
and 31 GHz and located its peak at ~ 0.8. Radio source properties 
in cluster fields have been characterized in e.g., Coble et al. (2007) 
tend to have a steep spectrum. In particular, the 1,4-to-31 GHz spec¬ 
tral index for the Coble et al. study had a mean value of 0.72. Sayers 
et al. (2013) explored the 1.4-to-31 GHz radio source spectral prop¬ 
erties towards 45 massive cluster systems and obtained a median 
value for a s of 0.89, which they showed was consistent with the 
30-to-140GHz spectral indices. The radio source population used 
to estimate the contamination to the Planck 143 GHz signal is likely 
to be a combination of field and cluster-bound radio sources due to 
the size of the Planck beam and the fact that some of the candidates 
might be spurious Planck detections. Overall, given the differences 
in the source selection and in frequency, and the agreement with 
other studies, our choice for of a s = 0.72 seems to be a reasonable 
one. 

In Table 10 we list the sum of all the predicted radio-source- 
flux densities at 100 and 143 GHz of all the NVSS-detected sources 
within 5' of our pointing centre. This yields an approximate mea¬ 
sure of the radio-source contamination in the Planck beam at these 
frequencies. The mean of the sum of all integrated source-flux 
densities at 1.4GHz is 61.0mJy (standard deviation, s.d. 71.6); 
at 100GHz it is 2.8 mjy (s.d.= 3.3) and at 143 GHz it is 2.2 mjy 
(s.d.=2.6). The SZ decrement towards each cluster candidate within 
the 143 GHz Planck beam is given in Table 10, together with 
the (expected) percentage of radio-source contamination to the 
Planck cluster signal at this frequency, which on average amounts 
to ~ 2.9%. The mean percentage contamination to the Planck SZ 
decrement would drop to ~ 1.3% if we used the Sayers et al. (2013) 
a s = 0.89 and would increase to ~ 5% if we used a flatter a s of 0.6. 
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Table 9. F 500 and 6500 values derived from the release 1 Planck all-sky maps . The second and third columns contain F 500 and 6500 from a blind analysis of the 
Planck maps, that is, from an analysis of Planck data alone, without any constraints from ancillary data. The fourth column contains the Planck F500 values 
calculated at the CARMA -8 cluster centroids using the CARMA- 8 -derived 6500 measurements. The ratio of F 500 values from columns 2 and 4 is given in the 
fifth column, while the ratio of 6500 from the blind Planck (column 2) and the CARMA -8 analyses is provided in column 6 (where the CARMA -8 6500 values 
have been taken from Table 3. 


Cluster ID 

F500 blind 

0500 blind 

T500 

l500,blmd/T500 

©5OO,blind/05OO 


arcmin 2 

arcmin 

arcmin 2 




x 10~ 4 


x 10~ 4 



POM 

13.2 ±6.4 

4.23 

8-2 

1.61 

0.98 

P086 

6.9 ± 2.5 

4.23 

5.8^;’ 

1.18 

1.14 

P097 

3.8 ±0.8 

0.92 

n 1 +0.6 

1.22 

0.29 

P109 

5.4 ± 1.2 

0.92 

5-9±}f 

0.91 

0.30 

P170 

11.7 ±4.1 

4.75 

7.1«;° 

1.65 

1.76 

P187 

11.5 ±4.1 

3.35 

1°-5±|;X 

1.10 

0.82 

P190 

7.9 ±4.7 

4.75 

5.9±i;» 

1.33 

1.36 

P205 

10.3 ±4.9 

3.35 

io-7±| ; | 

0.97 

0.68 

P351 

14.2 ±9.7 

6.75 

8.9j£? 

1.60 

1.27 


Table 10. Estimated radio-source contamination in Planck for clusters in our sample. Column 2 provides the sum of all the 1,4-GHz ‘deconvolved’ integrated 
flux-density measurements of NVSS sources detected within 5' of the Planck cluster centroid. Assuming a spectral index, a s , of 0.72 (see Paper 1 for details), 
we extrapolated the NVSS flux densities to find the total flux density at 100 and 143 GHz (columns 3 and 4), the most relevant Planck bands for SZ. The fifth 
column contains the (candidate) SZ decrement measured in the Planck 143-GHz band within its beam area in mJy. The last column presents the percentage 
radio-source contamination to the 143-GHz Planck candidate SZ decrement. This percentage is likely to be the amount by which the SZ flux is underestimated 
in the Planck data, since these are faint sources below the Planck point source detection limit. 


Cluster 

ID 

£ Flux 

densities 

at 1.4 GHz 

£ Flux 

densities 

at 100 GHz 

£ Flux 

densities 

at 143 GHz 

143-GHz Planck 

SZ decrement 

inside Beam 

% Radio Source 

Contamination 

to Planck SZ 

P014 

131.20 

6.10 

4.70 

- 77.3 

6.1 

P028 

120.70 

5.60 

4.40 

-90.2 

4.9 

P031 

14.10 

0.60 

0.50 

-57.6 

0.9 

P049 

62.60 

2.90 

2.30 

-55.7 

4.1 

P052 

293.00 

13.60 

10.50 

-67.6 

15.5 

P057 

126.80 

6.00 

4.60 

-96.5 

4.8 

P090 

4.30 

0.20 

0.20 

-42.0 

0.5 

P097 

5.90 

0.20 

0.20 

-38.4 

0.5 

P109 

54.90 

2.60 

2.00 

-108.4 

1.8 

P121 

31.30 

1.50 

1.10 

-79.9 

1.4 

PI 34 

15.80 

0.70 

0.60 

-54.9 

1.1 

P138 

37.00 

1.70 

1.30 

-66.8 

1.9 

P170 

18.40 

0.90 

0.70 

-141.1 

0.5 

P187 

71.10 

3.40 

2.60 

-70.4 

3.7 

P190 

18.90 

0.80 

0.70 

-48.6 

1.4 

P205 

15.80 

0.70 

0.60 

-82.1 

0.7 

P264 

8.40 

0.40 

0.30 

-61.1 

0.5 

P351 

67.10 

3.10 

2.40 

-101.7 

2.4 


Thus, we expect the flux density from unresolved radio sources 
towards our cluster candidates to be an insignificant contribution 
to the Planck SZ flux although individual clusters may have radio 
source contamination at the ~ 5 — 15% level. 


3.3 Improved Constraints on I 500 and 0 5 qo from the Use of a 
Planck Prior on Fsoo in the Analysis of CARMA -8 Data 

Due to their higher resolution (a factor of > 5), the CARMA-8 data 
are better suited than the Planck data to constrain 6500 - On the other 
hand, the large Planck beam (FWHMs; 10' at 100 GHz) allows the 


sampling parameter Yj for our clusters (all of which have 6)500 5') 

to be measured directly, which is not the case for the CARMA -8 
data due to its finite sampling of the uv plane and the missing zero¬ 
spacing information (a feature of all interferometers). We have ex¬ 
ploited this complementarity of the Planck and CARMA -8 data to 
reduce uncertainties in F 500 and 6)500 ■ I n order to do this, we fil¬ 
tered out the parameter chains (henceforth chains) for the analysis 
of the CARMA -8 data (model I) that had values of I 500 outside 
the range allowed by the Planck F 500 results (Table 9). We refer to 
the results from the remaining set of chains as the joint results (Ta¬ 
ble 13). In Figure 5 we plot the 2D marginalized distributions for 
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Table 11. Mean and 68% confidence uncertainties for radio-source parameters for sources detected in the LB CARMA-8 data towards candidate Planck 
clusters with a CARMA-8 SZ detection. These parameters have been obtained from the joint cluster+sources fit in McAdam to the SB data using cluster 
model I (Section 2.2.1). 


Source ID 

Cluster ID 

RAj/deg 

Dec 5 /deg 

S 3 i/Jy 


1 

POM 

240.831 ±0.001 

3.282 ± 0.001 

0.0081 ±0.0012 

0.5 ± 0.4 

2 

POM 

240.875 ± 0.001 

344.155± 0.001 

0.0089 ± 0.0022 

0.6 ± 0.4 

1 

P109 

275.718 ±0.002 

78.384 ± 0.002 

0.0018 ±0.0004 

0.6 ± 0.5 

1 

P170 

132.813 ±0.002 

48.619 ± 0.002 

0.0046 ± 0.0007 

0.5 ± 0.5 

1 

P187 

113.084 ±0.001 

31.688 ±0.001 

0.0037 ± 0.0004 

0.5 ± 0.5 

1 

P351 

226.077 ± 0.002 

-5.914 ± 0.002 

0.0028 ± 0.0005 

0.6 ± 0.5 


Table 12. Bayes Factor, K, for SZ signals detected in the CARMA-8 data 
by McAdam. Since the prior ratio is set to unity, the Bayes Factor provides 
a measure of the quality of the model fit to the data. For two clusters, POM 
and PI34, two potential SZ signals were detected in the Field of View (FoV) 
of each observation; the second cluster-like signal is labelled with a ‘b\ 
We adopt Jeffreys (1961) interpretation of K, though with fewer categories, 
to be even more conservative. We consider K ^ 0.1 to be strong evidence 
against a cluster signal (NC); 0.1 < K ^ 10 means our data cannot be used 
on their own to distinguish robustly between a model with or without a 
cluster signal (ND) and K > 10 indicates there is strong evidence for the 
presence of a cluster signal in the data (D). For reference, the signal-to- 
noise ratio (SNR) for detected clusters (those highlighted in bold font) in 
CARMA-8 and Planck data are given in Table 2 of Paper 1. 


Cluster 

Bayes Factor 

Degree of detection 

P014 

3.4e+01 

D 

POMb 

4.2e-01 

ND 

P028 

4.6e-01 

ND 

P031 

5.0e-02 

NC 

P049 

2.2e-01 

ND 

P052 

9.5e-01 

ND 

P057 

6. le-01 

ND 

P086 

3.5e+03 

D 

P090 

1.2e+00 

ND 

P097 

7.9e+02 

D 

P109 

8.4e+02 

D 

P121 

4.9e-01 

ND 

P134 

3.0e-02 

NC 

PI 34b 

3.0e-04 

NC 

P170 

7.2e+08 

D 

P187 

8.0e+08 

D 

P190 

5.6e+18 

D 

P205 

1.3e+09 

D 

P264 

2.3e+00 

ND 

P351 

2.7e+01 

D 


T 500 and 6500 for the CARMA -8 data alone (black contours) and for 
the joint results (magenta contours). Similar approaches comparing 
Planck data with higher resolution SZ data have been undertaken 
by Planck Collaboration et al. (2013e) (with AMI). Muchovej et 
al. (2012) (with CARMA), Sayers et al. (2013) (with BOLOCAM) 
and AMI and Planck Consortia: Perrott et al. (2014) (with AMI). 
Clearly, the introduction of cluster size constraints from high res¬ 
olution interferometry data provides a powerful way to shrink the 
uncertainties in I 500 — 6500 phase space. 


4 DISCUSSION 
4.1 Use of Priors 

When undertaking a Bayesian analysis, it is important not only 
to check that the priors on individual parameters are sufficiently 
wide, such that the distributions are not being truncated, but also 


Table 13. Columns 2 and 3 contains the joint 6500 and F 500 values, which 
were computed by truncating the output of the CARMA -8 chains to have 
T 500 values in the range allowed by the Planck results (column 4 in Table 

9). 


Cluster ID 

®500, Joint 

^500,Joint 


arcmin 

arcmin 2 



x 10 -4 

POM 

3-8±8i 

8-2±g;« 

P086 

7 7 +O .3 
z,/ -0.4 

5 . 518:1 

P097 

2 1+ 0 - 1 
z,a -0.3 

•7 -7+0.2 
-0.2 

P109 

7 7 +O .2 
^ -0.4 

5 7 + 0 - 3 
J,/ -0.5 

P170 

•7 -7+0.2 

7 l +0 - 4 

P187 

4 7+0-3 
4 - z -0.4 

1Q+ 10 

1U — 1.4 

P190 

7 o+0.2 
Z -°-0.2 

C O+0.5 

^■o-o. 7 

P205 

4 3 +0-3 
^-o.s 

10+ 11 

1 —1.5 

P351 

C c+ 0.2 

0.8 

0 +O . 7 
— 1.2 


that the effective prior is not biasing the cluster parameter results. 
Here the term effective prior refers to the prior that is being placed 
on a model parameter while taking into account the combined ef¬ 
fect from all the priors given to the set of sampling parameters. 
What may seem to be inconspicuous priors on individual parame¬ 
ters can occasionally jointly re-shape the high dimensional param¬ 
eter space in unphysical ways; this was noticed in e.g., AMI Con¬ 
sortium: Zwart et al. (2011). Biases from effective priors should be 
investigated by undertaking the analysis without data i.e., by set¬ 
ting the likelihood function to a constant value. Such studies for the 
models used in this work have been presented in AMI Consortium: 
Rodrfguez-Gonzalvez et al. (2012), AMI Consortium: Olamaie et 
al. (2012) and Olamaie, Hobson, & Grainge (2013) and have deter¬ 
mined that the combination of all the model priors does not bias the 
results. 


4.2 Characterization of the Cluster Candidates 

4.2.1 Cluster position and Morphology 

The mean separation (and standard deviation, s.d.) of the CARMA- 
8 centroids from Model I and the Planck position is ~ lA' (0.5); see 
Tables 3 and 6 for offsets from the CARMA -8 SZ decrement to the 
Planck position. This offset is comparable to the offsets between 
PLANCK and X-ray cluster centroids found for the ESZ (Planck 
Collaboration et al. 2011 VII) and the PSZ (Planck Collaboration 
et al. 2013 XXIX), which were typically 2' and 70", respectively. 
The cluster candidate with the largest separation, ~ 2.5', is P014. 
The high-resolution CARMA -8 data allow for the reduction of po- 
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Figure 2. Average 6500 . ( 05 oo}> values calculated from the MCXC cata¬ 
log (Piffaretti et al. 2011) within a a series of redshift bins, as indicated 
in the x-axis e.g., the first bin represents (6500) f° r MCXC clusters with 
0.1 ^ z < 0.2. The barplots filled with grey diagonal lines use all the rel¬ 
evant MCXC cluster entries to get (6500). while those filled with blue di¬ 
agonal lines, only include the more massive clusters, M 500 > 4 x 10 14 Mq, 
which we expect to be more representative of our target sample as Planck 
detects the most massive clusters. Standard deviations are displayed as ver¬ 
tical lines. The average CARMA- 8 -derived 6500 estimate for our clusters is 
depicted with an orange horizontal line. This plot suggests the typical 6500 
values for our clusters are most comparable with the ±00 values for MCXC 
clusters at 0.2 < z < 0.4. Furthermore, when combined with our photomet¬ 
ric redshift measurements in Paper 1 which place our clusters at z ~ 0.5, 
the data suggest that we are finding the largest, most-massive, clusters at 
intermediate redshifts. 


sitional uncertainties in the Planck catalog for candidate clusters 
from a few arcminutes to within < 30”. This is crucial, amongst 
other things, for the efficient follow-up of these candidate systems 
at other wavelengths. P351 has the largest positional uncertainties 
for both parameterizations (~ 40”), indicative of a poorer fit of the 
models to the data, since the noise in the CARMA-8 data is one 
of the smallest of the sample. Interestingly, this cluster stands out 
in the /3 parameterization for having the shallowest profile (Figure 
1), and in the gNFW parameterization for having the largest 6500 - 
Overall, the positional uncertainties from the shape-fitting model I, 
tends to be larger than that from the radial profile based model II, 
typically by a factor of 1.2 and reaching a factor of 2.8. The differ¬ 
ent parameter degeneracies resulting from each analysis is likely 
to be the dominant cause for this. As shown in Figure 1 of AMI 
Consortium: Rodrfguez-Gonzalvez et al. (2012), in the Tgoo — 6500 
plane, the 2D marginalized distribution for the cluster size is sig¬ 
nificantly narrower for the ft parameterization (model II) than for 
the gNFW parameterization (model I). 

On average, cluster candidates with CARMA -8 detections 
have Q s = 3.4±0.2' with an s.d. of 0.9 and 6500 = 3.9± 1.6\ with 
s.d. of 0.9 (see Table 3). The largest cluster has 6 S = 5', 6500 = 5.3 r 
(P351) and the smallest 6 S = 2', 0gQ () = 2.1' (P170). In Paper 1 we 
estimated the photometric redshifts for our cluster candidates with a 
CARMA -8 SZ detection and found that, on average, they appeared 


to be at z ~ 0.5; in Section 4.3.1 we report on the spectroscopic 
confirmation of P097 at z = 0.565. The relatively small values for 
0500 would support the notion that our systems are at intermediate 
redshifts (z > 0.5). In comparison, for the MCXC catalog of X-ray- 
identified clusters (Piffaretti et al. 2011), whose mean redshift is 
0.18, the mean X-ray-derived 6500 is a factor of 2 larger. In Figure 
2 we plot the average @500 within a series of redshift ranges start¬ 
ing from z = 0.1 for all MCXC clusters (in grey) and for only the 
more massive, M 500 > 4 x 10 14 Mq, clusters (in blue), which should 
be more representative of the cluster candidates analyzed here (see 
Figure 3) and mark the average CARMA- 8 -derived 0goo for our 
clusters with an orange line. This plot suggests the ©500 values for 
our clusters are most comparable with the @500 values for MCXC 
clusters at 0.2 < z < 0.4. 

The resolution of the CARMA-8 data, together with the often 
poor signal-to-noise ratios and complications in the analysis, e.g., 
regarding the presence radio sources towards some systems, makes 
getting accurate measurements of the ellipticity r] challenging, with 
typical uncertainties in p of 0.2 (Table 3). The mean and standard 
deviation of 77 for our sample is 0.7 ±0.07. The values are there¬ 
fore consistent with unity, to within the uncertainties. Nevertheless, 
the use of a spherical model is physically motivated and allows the 
propagation of realistic sources of uncertainty. Moreover, compar¬ 
ison of models with spherical and elliptical geometries for similar 
data from AMI is presented in AMI Consortium: Hurley-Walker 
et al. (2012) which show the Bayesian evidences are too alike for 
model comparison, indicating that the addition of complexity to the 
model by introducing an ellipticity parameter is not significantly 
penalised. In AMI and Planck Consortia: Perrott et al. (2014) mod¬ 
elling of AMI cluster data with and ellipsoidal GNFW profile in¬ 
stead of a spherical profile had a negligible effect on the constraints 
in Y. Our CARMA data with higher noise levels and generally more 
benign source environments should show even smaller effects. 

T] values close to 1 would be expected for relaxed systems, 
whose projected signal is close to spheroidal, unless the main 
merger axis is along the line of sight. On the other hand, disturbed 
clusters should have rj —> 0.5. Some evidence for a correlation of 
cluster ellipticity and dynamical state has been found in simula¬ 
tions, e.g, Krause et al. (2012) and data, e.g., Kolokotronis et al. 
(2001) (X-ray), Plionis (2002) (X-ray and optical) and AMI Con¬ 
sortium: Rodrfguez-Gonzalvez et al. (2012) (SZ), although this cor¬ 
relation has a large scatter. Hence, from the derived fits to the data, 
we conclude that our sample is likely to be mostly comprised by 
large, dynamically active systems, unlikely to have fully virialized 
which is not surprising given the intermediate redshifts of the sam¬ 
ple. 

4.3 Cluster-Mass Estimate 

To estimate the total cluster mass Mgoo within 7 - 500 , we use the Ola- 
maie, Hobson, & Grainge (2013) cluster parameterization, which 
samples directly from Moqq 6 . This model describes the cluster dark 
matter halo with an NFW profile (Navarro, Frenk, & White 1996) 
and the pressure profile with a gNFW profile (Nagai, Kravtsov, & 

6 M 500 i s determined by calculating r$ 00 , which in turn is computed by 
equating the expression for mass from the NFW density profile within r^m 
and the mass within an spherical volume of radius rsoo under the assumption 
of spherical geometry. Under this NFW-gNFW-based cluster parameteriza¬ 
tion the relation between M 200 and M 50 O is M 200 = 1.35 x M 5 oo- For further 
information on this cluster model see Olamaie, Hobson, & Grainge (2013). 
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Figure 3. M500 estimates for candidate cluster P190 as a function of red- 
shift. To obtain the M500 measurements a third parameterization, model III 
(Section 4.3; Olamaie, Hobson. & Grainge 2013). which samples directly 
from the cluster mass, was implemented in the McAdam analysis. Since this 
model requires redshift information as an input, but spectroscopic-redshift 
information is not available for most of our candidate clusters, we ran this 
model run six times, from z = 0.1 to z — 1.0 in steps of 0.2. The range of 
M500 values encompassing 68% of the probability distribution are shown in 
the red shaded area. We overplot the M500 estimates from the Planck PSZ 
cluster catalog (Planck Collaboration et al. 2013 XXIX) for entries with an 
associated redshift in blue + signs. The PSZ catalog of clusters contains the 
majority of the most massive X-ray-detected clusters in the MCXC cata¬ 
log. Hence, the clusters in our sample, which have SZ signatures similar to 
that of PI90, are amongst the most massive known systems at intermediate 
redshifts. 



Figure 4. Keck/MOSFIRE Y-band image for P097 in grey background, 
where high intensity regions appear darker. Overlaid are the CARMA SZ 
contours, as shown in Paper I, in blue. The contours are 2, 3 and 4(7, 
where (7, the RMS of the short baseline CARMA data towards P097, is 
0.653 mJy/beam. The red circle shows the galaxy the spectroscopic target 
whose redshift was determined to be 0.565 


Vikhlinin 2007), using the set of gNFW parameters derived by Ar- 
naud et al. (2010). There are two other additional assumptions of 
this model: ( 1 ) the gas is in hydrostatic equilibrium and ( 2 ) the 
gas mass fraction is small compared to unity (< 0.15). 7 The clus¬ 
ter redshift is a necessary input to this parameterization. In the ab¬ 
sence of spectroscopic data towards our cluster candidates to get an 
accurate redshift estimate 8 * * , coarse photometric redshifts based on 
SDSS and WISE colors were calculated in Paper 1. In Figure 3 we 
plot the M 500 estimate (mean values are depicted by the thick line 
and the area covering the 68 % of the probability distribution has 
been shaded in red) as a function of z for one of our cluster can¬ 
didates, PI90. To produce this plot we ran the Olamaie, Hobson, 
& Grainge (2013) cluster parameterization six times using a Delta- 
prior on redshift, which we set to values from 0.1 to 1.0 in steps 
of 0.2. We chose P190 since it is quite representative of our cluster 
candidates and, at 8(7, where (7 is the short baseline rms, it has the 
best SNR of the sample (see Table 2 in Paper 1). The photometric- 
redshift estimate for P190 from Paper 1 was 0.5, which is also the 
average expected photometric redshift for the sample of CARMA -8 
SZ detections. At this redshift, M 500 = 0.8±0.2 x 10 15 M@. As seen 
in Figure 3, after z ~ 0.3 M 500 is a fairly flat function of z (since 
the z dependence in the model is carried by the angular diameter 
distance), such that, to one significant figure, our mean value is 
identical. 


4.3.1 Spectroscopic Redshift Determination for P097 

We have measured the spectroscopic redshift of a likely galaxy 
member of P097 through Keck/MOSFIRE Y-band spectroscopy. 
We deem it a likely member, given that it is situated close to the 
peak of the CARMA SZ decrement (within the 4 0 contour) and 
close to a group of tightly clustered galaxies, as shown in Figure 4, 
of similar colours. We detect clear evidence for Ha, [Nil], [SII] at 
redshifts of 0.565, with the strongest lines being present in a Sloan- 
detected galaxy located at RA,DEC of (14:55:25.3,+58:52:33.86; 
J2000). We also find evidence of velocity structure in this galaxy 
with the Ha being double peaked with the line components sep¬ 
arated by ~150 km/s. These data will be published in a separate 
paper. We calculated M 500 for this cluster, as we did for PI90 in 
the previous section, setting the redshift prior to a delta function 
at z = 0.565 and obtained M 500 = 0.7 =1=0.2 x 1O 15 M 0 , supporting 
the notion that our sample of clusters are some of the most massive 
clusters at z > 0.5. Further follow-up of this sample in the X-rays 
and through weak lensing measurements with Euclid will help con¬ 
strain the mass of these clusters more strongly. 

4.4 The F 5 OO- 05 OO Degeneracy 

In Figure 5, the 2D marginalized posterior distributions in the 
F 500 — 0500 plane for (i) the CARMA -8 data alone and (ii) the joint 
analysis of CARMA -8 and Planck data are displayed for each clus¬ 
ter candidate with an SZ detection in the CARMA -8 data. It is clear 

7 Extensive work has been done to characterize the hydrostatic mass bias, 
which has been shown to range between ~ 5 and 45% depending on dy¬ 
namical state in numerical simulations (e.g., Rasia et al. 2006, Nagai et al. 
2007b and Molnar et al. 2010) and studies comparing SZ and weak lensing 
mass estimates (e.g., Zhang et al. 2010 and Mahdavi et al. 2013). 

8 PI87 might an exception as the observational evidence shows that it is 

likely to be associated with a well-known Abell cluster, Abell 586, at z = 

0.171. We obtained spectroscopic confirmation for P097, see Section 4.3.1 
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(a) Left to right: POM, P086, P097 



(b) Left to right: P109, P170, P187 





(c) Left to right: P190, P205, P351 


Figure 5. 2D posterior distributions in the F 500 — $500 plane for the CARMA-8 data alone in black contours and for the CARMA-8 data using a Planck prior 
on Y 500 (column 4 in Table 9), in magenta contours, to which we refer to as the joint CARMA -Planck constraints. The y-axis is Y 500 in units of arcmin 2 and 
the x-axis is $500 in units of arcmin. The lower and upper limits of the y-axis are 0 and 0.006 arcmin 2 , labelled in steps of 0.001 arcmin 2 and for the x-axis 
they are 0 to 12 arcmin, labelled in steps of 2 arcmin. The inner and outer contours in each set indicate the areas enclosing the 68% and 95% of the probability 
distribution. For all the candidate clusters there is dramatic improvement in F 500 — O 500 uncertainties when the CARMA and Planck data are jointly analyzed 
to yield cluster parameters. 
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Figure 6. Left: F500 measured by CARMA -8 against the blind Planck No strong bias indications are detected between the Planck, blind and CARMA -8 
F500 results. The best-fit (red, dotted) and 1:1 (black, solid) lines are included in this plot. Right: ratio of the Planck, blind F500 and the joint F500 values against 
the ratio of the CARMA -8 F500 and the joint Fsqo- For F500, both the CARMA -8 and Planck, blind values are practically always higher than the joint values by 
as much as a factor of 2 . 
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Figure 7. Left: @500 measured by CARMA-8 against the blind Planck $500 • Uncertainties on the Planck, blind #500 are not included as they are determined by 
the large FWHM of the Planck beam and can range up to 10 '. The best-fit (red, dotted) and 1:1 (black, solid) lines are included in this plot. Right: Ratio of the 
blind Planck 6500 and the joint (Ton against the ratio of the CARMA-8 (Too and the joint (Too- No strong systematic differences are seen between the different 
6500 estimates. The higher resolution of the CARMA -8 data allows ©500 to be constrained much more strongly than using Planck data alone. 


from these plots and from Table 9 that there is good overlap be¬ 
tween the Planck and CARMA -8 derived cluster parameter space 
for all of the candidate clusters. The range of Planck T 500 values 
after application of the 6500 priors from CARMA -8 are within the 
68 % contours for the CARMA-only analysis. While the Planck 


only range of T 500 ,Mind values is wide, in combination with the 
0500 constraints from the higher resolution CARMA-8 data, the 
I 500 — 0500 space is significantly reduced. To explore this further, 
in Figure 6 we have plotted: 

(i) the relation between T 500 from CARMA-8 and from the Planck, 
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Table 14. This table shows the average improvement in the SNR of the 
I 500 measurements taken from the independent analysis of (blind) Planck 
or C ARM A-8 data versus the joint analysis of Planck and C ARM A-8 
data. For clarity in the table we use a shorthand notation. We refer to 
SNRy 500 Joint as the ratio of F500,Joint and the 68% error in F 50 o, Joint, that is, 
I500,joint/< tF>oo,J oint, where, as denoted by the suffix label, the results are ob¬ 
tained from the joint analysis of the CARMA-8 and Planck data. Similarly, 

SNR J'500,bM 311(1 SNR l500,CARMA-8 m ' e ^ ratlOS ° f F 500 and tlle 68% elTOT 

in Y 500 for the independent analysis of Planck (blind) data and CARMA- 
8 data, respectively. Since we have asymmetric errors for F 500 , in the first 
column we provide the SNR ratios calculated using the lower-bound errors 
and in the second column those calculated using the upper-bound errors. In 
an entirely analogous manner, we calculate the improvements in the SNR 
of 6500 ns we shift from independent analyses of the Planck or CARMA -8 
data to a joint analysis. In the lower sub-table, we focus on the improve¬ 
ments to the 68%-confidence uncertainties in I 500 and (boo from undertak¬ 
ing joint Planck and CARMA-8 analyses rather than using each data set on 
its own. It is clear that a joint analysis of the two datasets provides much 
tighter constraints in the I^oo — 6500 parameter space. This is mostly due to 
the challenge in constraining cluster sizes in the Planck data and a lack of 
sensitivity to the signal from the outskirts of the cluster in the CARMA-8 
data. 



trAoQ.CARMA- 8 ' 

C ^ 800 ,Joint I 

g %00,CARMA-8 

ff0 5OO, Joint 
/ g f5QQ,blind \ 

\ 1,^500,Joint / 

’ CT0 5OO,Joint ' 


6.9 9.4 
7.3 4.4 
7.3 5.2 
52.2 27.7 


blind analysis Fgoo, blind (left panel) 

(ii) the ratio of Tgoo from the Planck, blind analysis and from the 
joint results against the ratio of Y 50 O from fits to CARMA -8 data 
and from the joint analysis (right panel). 

In Figure 7, we present similar plots for 6500 - Inspection of Fig¬ 
ures 6 and 7, left panels, shows that the Planck, blind values for 
K 500 , as well as for 6500 , appear unbiased with respect to those 
from CARMA- 8 . As expected, due to Planck’ s low angular reso¬ 
lution, the overall agreement between the Planck, blind measure¬ 
ments and those of CARMA -8 are much better for I 500 than for 
0500. wit h (©500,CARMA-8/©500, Planck blind) = 1-5, S.d.= l.l and 
( f 500,CARMA-8/ F 500, Planck blind) = 1-1- s.d.=0.4. This good agree¬ 
ment in I 500 is in contrast with recent results by Von der Linden et 
al. (2014), who compare cluster-mass estimates derived by Planck 
and weak lensing data towards 22 clusters and find Planck masses 
are lower than the weak lensing masses typically by ~ 30%. Such 
a bias would indeed alleviate the tension found in (Planck Collab¬ 
oration et al. 2013 XX), where the 95% probability contours in the 
tfg — £2 m plane derived from cluster data and from the CMB tem¬ 
perature power spectrum do not agree, when accounting for up to 
a 20 % bias from the assumption of hydrostatic equilibrium in the 
X-ray-based cluster-mass scaling relations. Since we find no signs 
of bias between the I 500 measurements from Planck and C ARMA- 
8 , this might indicate the bias arises when comparing masses rather 
than I 500 , that is, from the choice of scaling relations used to es¬ 
timate the cluster mass. Although, our cluster sample is relatively 
small and our parameter uncertainties are substantial other sources 
of bias in the SZ, X-ray and lensing measurements need to be inves¬ 
tigated further with the same samples of objects rather than cluster 
samples selected in different ways and extending over different red- 
shift ranges. 

The right panel of Figure 6 shows that, for all but two candi¬ 
date clusters, the K 500 measurements derived from the joint analy¬ 
sis are consistently lower than from the independent analysis of the 
Planck, blind and CARMA -8 datasets, sometimes by as much as 
a factor of 2. The average magnitude and range of the shift from 
the independent-to-joint I 500 values are similar for the CARMA- 
8 and the Planck, blind analysis, with / y 0Q,1>1 ‘ nd \ = 1.3, s.d.= 0.3, 

( y ToZ,t~ 8 ) = L4 ’ s - d - = °- 4 ' In the case of 0 5OO, Figure 7 (right), 
there is no systematic offset between the joint results and those 
from either CARMA -8 or Planck, blind. On average, the agree¬ 
ment between the joint and independent results appears to be good, 
/gSOflaA = ! 0 s d _ 0.5 / 0 5OO,CARMA-8 \ = J J S ( J = Q.2, but, 

\ ^(H),Joint / \ '^'500,Joint / 

in the case of the Planck, blind measurements, the large standard 
deviations are an indication of its poor resolution. 

In Table 14 we quantify the improvement in the constraints for 
F 500 and 6)500 derived from the independent analyses of Planck and 
CARMA -8 data with respect to the joint results. The largest im¬ 
provement is seen for 6)500 Planck blind, as this parameter is only 
weakly constrained by Planck data alone, with an associated uncer¬ 
tainty anywhere up to ~ 10'. Yet, the advantage of a joint Planck 
and CARMA -8 analysis is very significant both for 6500 and F 500 
and for both the CARMA -8 and Planck results, with uncertainties 
dropping by > 400%. As mentioned in Section 3.3, the SZ mea¬ 
surements from Planck and CARMA -8 data are complementary as 
they probe different cluster scales at different resolutions. More¬ 
over, since the F 500 — 6)500 degeneracy for each data set have differ¬ 
ent orientations (an effect already reported in Planck Collaboration 
et al. 2013e), a joint analysis looking at the overlapping regions will 
result in a further reduction of parameter space. 
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Figure 8. Left: Plot of I 500 derived from two different analyses of PLANCK data against I 500 derived from CARMA -8 data. We plot two sets of >500 based on 
Planck measurements (i) using Planck data alone; these data points are referred to as ‘blind’ and are shown in red crosses (ii) using the range of F 500 Planck 
(blind) values to constrain further the *500 values obtained in the analysis of CARMA -8 data; these are referred to as 'joint' values and are shown in green 
crosses inside a square. In addition, we plot a 1 : 1 line in solid black. For the joint F 500 values, these are plotted on the y-axis, whilst keeping the x-axis 
as the CARMA-8 I 500 values. Right: Plot of F 500 from several analyses of PLANCK data against F 500 from AMI. We plot two sets of F 500 based on Planck 
measurements (i) using Planck data alone; these data points are referred to as 'blind' and are shown in blue diamonds (ii) using the range of F 500 Planck 
(blind) values as a prior in the analysis of the AMI data; these are referred to as 'joint' values and are shown in purple + signs inside a diamond. For the joint 
*500 values, these are plotted on the y-axis, whilst keeping the x-axis as the AMI F 500 values. In addition, we plot a 1 : 1 line in solid black. Comparison of 
the left and right plots shows the Planck (blind) measurements are generally in good agreement with CARMA-8, with no signs of systematic offsets between 
the two measurements. Calculating the Planck F 500 using the CARMA -8 6500 and position measurements, decreases the Planck F 50 O estimates, which become 
systematically slightly lower than the CARMA -8 estimates (Table 9). This systematic difference is enhanced for the joint F 500 values. For AMI, the Planck 
blind F 500 values appear to be consistently higher than the AMI values for all but one cluster. While this difference narrows for the joint AMl-Planck F 500 
estimates, it is not resolved. Hence, the AMI F 500 estimates are generally higher than those for Planck, irrespective of the choice of Planck -derived F 500 while 
CARMA -8 F 500 values are in good agreement with the (blind) Plcmck results, yet are generally higher when priors are applied to the CARMA-8 or Planck 
data. 


4.5 Comparison with the AMI -Planck Study 

The Arcminute MicroKelvin Imager (AMI Consortium: Zwart et 
al. 2008) has followed up ss 100 Planck -detected systems, most 
of which are previously confirmed systems. Comparison of AMI 
and Planck results for 11 clusters in Planck Collaboration et al. 
2013e showed that, as seen by Planck, clusters appear larger and 
brighter than by AMI. This result has now been confirmed in a 
larger upcoming study (AMI and Planck Consortia: Perrott et al. 
2014). In Figure 8 we compare AMI and CARMA -8 values for Fsoo 
against the PLANCK results. We note that the analysis pipeline for 
the processing of CARMA -8 data in this work is identical to that of 
AMI, allowing for a clean comparison between both studies. 

All eleven clusters in the AMI-Planck study are known X- 
ray clusters (except for two) and are at 0.11 < z < 0.55, with an 
average z of 0.3 and a mean 0 5 oo of 4.8. Our sample of cluster can¬ 
didates is expected to have a mean (coarse) photometric redshift of 
~ 0.5 (see Paper 1) and a mean 6500 of 3.9. The smaller angular 
extent of our sample of objects could be indicative that they are in 
fact at higher redshifts. As pointed out in Planck Collaboration et 
al. (2013e), the F 500 measurements from Planck are systematically 


higher than those for AMI, (Fsoo,AMl/*500,blindPlanck) = 0.6, s.d. 
= 0.3. For the CARMA -8 results, on the other hand, we find good 
agreement between the CARMA -8 and Planck -blind F 500 values, 
( *500,CARMA—8 7*500,blind Planck) = 1-1. s -d- = 0.4, with no system¬ 
atic offset between the two measurements. 

Planck Collaboration et al. (2013e) point to the use of a fixed 
gNFW profile as a likely source for systematic discrepancies be¬ 
tween the AMI -Planck results and they plan to investigate changes 
to the results when a wider range of profiles are allowed in the fit¬ 
ting process. With relatively similar spatial coverage between the 
CARMA -8 and AMI interferometers, the impact of using a gNFW 
profile in the analysis of either data set should be comparable (for 
most cluster observations), though this will be investigated in detail 
in future work. The fact that, despite this, the CARMA -8 measure¬ 
ments are in good agreement with Planck could mean that either the 
CARMA -8 and Planck data are both biased-high due to systematics 
yet to be determined, or the AMI data are biased low, or a combina¬ 
tion of both of these options. We stress that the analysis pipeline for 
deriving cluster parameters in this work and in Planck Collabora¬ 
tion et al. (2013e) is the same. Apart from the (typically) relatively 
small changes in the uv sampling of both instruments, the major 








18 RocLriguez-Gonzdlvez et al. 


difference between the two datasets is that radio-source contamina¬ 
tion at 16 GHz tends to be much stronger than at 31 GHz, since ra¬ 
dio sources in this frequency range tend to have steep falling spec¬ 
tra and, hence, it tends to have a smaller impact on the CARMA -8 
data. That said, AMI has a separate array of antennas designed to 
make a simultaneous, sensitive, high resolution map of the radio¬ 
source environment towards the cluster, in order to be able to detect 
and accurately model contaminating radio sources, even those with 
small flux densities (~ 350/xJy/beam). 

A more likely alternative is the intrinsic heterogeneity in the 
different cluster samples. At low redshifts, when a cluster is com¬ 
pact relative to the AMI beam, heterogeneities in the cluster profiles 
get averaged out and the cluster integrated K 500 values agree with 
those from Planck. However, if the cluster is spatially extended rel¬ 
ative to the beam, a fraction of the SZ flux is missed and that results 
in an underestimate of the AMI SZ flux relative to the Planck flux. 
Since our sample is at higher redshifts and appears to be less ex¬ 
tended compared to the clusters presented in Planck Collaboration 
et al. (2013e), it is likely that profile heterogeneities are averaged 
out alleviating this effect seen in the low-z clusters. The prediction 
therefore, is that more distant, compact clusters that are followed 
up by AMI will show better agreement in 1500 values with Planck 
although hints of a size-dependency to the agreement are already 
seen in Planck Collaboration et al. (2013e). 

The results from this comparison between Planck, AMI and 
CARMA -8 SZ measurements draw further attention to the need to 
understand the nature of systematics in the data, in order to use ac¬ 
curate cluster-mass estimates for cosmological studies. To address 
this, we plan to analyze a sample of clusters observed by all three 
instruments in the future. 


5 CONCLUSIONS 

We have undertaken high (1-2') spatial resolution 31 GHz observa¬ 
tions with CARMA -8 of 19 Planck -discovered cluster candidates 
associated with significant overdensities of galaxies in the WISE 
early data release (> 1 galaxies/arcmin 2 ). The data reduction, clus¬ 
ter validation and photometric-redshift estimation were presented 
in a previous article (Paper 1; Rodrfguez-Gonzalvez et al. 2015). In 
this work we used a Bayesian-analysis software package to analyze 
the CARMA -8 data. Firstly, we used the Bayesian evidence to com¬ 
pare models with and without a cluster SZ signal in the CARMA -8 
data to determine that nine clusters are robust SZ detections and 
that two candidate clusters are most likely spurious. The data qual¬ 
ity for the remaining targets was insufficient to confirm or rule out 
the presence of a cluster signal in the data. 

Secondly, we analyzed the 9 CARMA -8 SZ detections with 
two cluster parameterizations. The first was based on a fixed-shape 
gNFW profile, following the model used in the analysis of Planck 
data (e.g., Planck Collaboration et al. 2011 VII), to facilitate a com¬ 
parative study. The second was based on a /3 gas density profile that 
allows for the shape parameters to be fit. There is reasonable corre¬ 
spondence for the cluster characteristics derived from either param¬ 
eterization, though there are some exceptions. In particular, we find 
that the volume-integrated brightness temperature within 6500 cal¬ 
culated using results from the /3 profile does not correlate well with 
T 500 from the gNFW parameterization for two systems. This sug¬ 
gests that differences in the adopted profile can have a significant 
impact on the cluster-parameter constraints derived from CARMA- 
8 data. Indeed, radial brightness temperature profiles for individual 
clusters obtained using the ft model results exhibit a level of het¬ 


erogeneity distinguishable outside the uncertainty, with the degree 
of concentration for the profiles within 6500 and within 600"(~ the 
FWHM of the 100-GHz Planck beam) being between ~ 1 and 2.5 
times different. 

Cluster-parameter constraints from the gNFW parameteriza¬ 
tion showed that, on average, the CARMA -8 SZ centroid is dis¬ 
placed from that of Planck by ~ 1.5 / . Overall, we find that our 
systems have relatively small 6)500 estimates, with a mean value 
of 3.9'. This is a factor of two smaller than for the MCXC clus¬ 
ters, whose mean redshift is 0.18. This provides further, tentative, 
evidence to support the photometric-redshift estimates from Paper 
1, which expected the sample to have a mean z of ~ 0.5. Using 
Keck/MOSFIRE Y-band spectroscopy, we were able to confirm the 
redshift of a likely galaxy-member of one of our cluster candidates, 
P097, to be at z = 0.565. 

We analyzed data towards PI90, a representative candidate 
cluster for the sample, with a cluster parameterization that samples 
from the cluster mass. This parameterization requires z as an input. 
We ran it from a z of 0.1 to 1 in steps of 0.2. Beyond the z > 0.3 
regime, the dependence of mass on z is very mild (with the mass 
remaining unchanged to within 1 significant figure). Our estimate 
for M 500 at the expected average z of our sample, 0.5 (Paper 1), is 
0.8 ±0.2 x 10 15 Mq . 

We compared the Planck (blind) and CARMA- 
8 measurements for I 500 and 6500 - Both sets of re¬ 
sults appear to be unbiased and in excellent agreement, 
with (*500,CARMA—8 /*500, blind,Planck ) = 1-1- S.d.=0.4. and 

(©500,CARMA-8/05OO, blind,Planck) = 1-5, S.d.=l.l, whose larger 
difference is a result of the poor spatial resolution of Planck. This 
is in contrast with the results from a similar study between AMI 
and Planck that reported systematic differences between these 
parameters for the two instruments, with Planck characterizing the 
clusters as larger (typically by rs 20%) and brighter (by ~ 35% 
on average) than AMI. However, it should be emphasized that the 
clusters studied in this paper are on average at higher redshift and 
more compact than those in the AMl-Planck joint analysis. 

Our results, and those from AMI, seem to not support the 
recent results by Von der Linden et al. (2014), who find Planck 
masses to be biased-low by ~ 30% with respect to weak-lensing 
masses and, potentially, with Planck Collaboration et al. (2013 XX) 
who consider the possibility of Planck masses to be biased-low by 
20 — 40% to explain inconsistencies in their results from cluster 
data and the CMB temperature power spectrum. The good agree¬ 
ment between Planck and CARMA -8 Y 500 measurements, in con¬ 
trast with the large differences in the Planck masses derived from 
X-ray-based scaling relations, could be an indication that the ori¬ 
gin of the discrepancy lies primarily in the choice of scalings and in 
the heterogeneity in cluster profiles with increasing redshift. How¬ 
ever, our sample size is small and our uncertainties substantial. A 
large multi-frequency study of clusters, which included Planck and 
a high-resolution SZ experiment, like CARMA -8 or AMI, to check 
for systematics in I 500 , as we ll as lensing and X-ray data to inves¬ 
tigate differences in cluster-mass estimates, would be beneficial to 
fully address this. 

We exploited the complementarity of the Planck and 
CARMA -8 datasets—the former measures the entire cluster flux 
directly unlike the latter, which can, on the other hand, constrain 
6 ) 500 — to reduce the size of the T 500 — 6)500 degeneracy by apply¬ 
ing a Planck prior on the T^oo obtained from CARMA -8 alone. We 
show how this joint analysis reduces uncertainties in *500 derived 
from the Planck and CARMA -8 data individually by more than a 
factor of > 5. 
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In this article and in its companion paper (Rodrlguez- 
Gonzalvez et al. 2015), we have demonstrated (1) an interest¬ 
ing technique for the selection of massive clusters at intermedi¬ 
ate z > 0.5 redshifts by cross-correlating Planck data with WISE 
and other data (2) a powerful method for breaking degeneracies in 
the > 500 (flux) — @ 5 oo(size) plane and thus greatly improving con¬ 
straints in these parameters. 
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